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Abstract 

The effect of coupling non-compact U(l) gauge fields to four dimensional 
simplicial quantum gravity is studied using strong coupling expansions and 
Monte Carlo simulations. For one gauge field the back-reaction of the matter 
on the geometry is weak. This changes, however, as more matter fields are 
introduced. For more than two gauge fields the degeneracy of random manifolds 
into branched polymers does not occur, and the branched polymer phase seems 
to be replaced by a new phase with a negative string susceptibility exponent 7 
and fractal dimension dn ~ 4. 



1 Introduction and motivations 

The statistical mechanics of random manifolds provides a possible framework for a 
non-perturbative construction of a quantum theory of gravity (there exist several 
excellent reviews, e.g. Q, [Q]). In this construction the standard recipes of lattice 
field theory are adopted, with a notable extension: the lattice itself becomes a dy- 
namical object instead of being an inert scaffold. For fixed topology, the summation 
over geometries involved in the partition function, approximating Feynman's path 
integral, is best implemented using the method of dynamical triangulations. The 
successes of this approach are particularly spectacular in two dimensions. The re- 
sults obtained with models exactly solvable in the continuum formalism have been 
reproduced. Furthermore, completely new results, hardly attainable with another 
approach, have been derived. As an example, let us quote the beautiful calculation 
of the "two-point" invariant correlation function in Ref. Q. The method of dynam- 
ical triangulations is also applicable to four dimensional gravity. However, and this 
is not a surprise, the construction of a viable theory is much more difficult in four 
than in two dimensions. 

At d = 2 the Einstein term in the action is a topological invariant and the dy- 
namics is insensitive to the corresponding coupling constant, as long as the topology 
of the manifold is kept fixed. For d > 2 one finds that changing the (lattice analogue) 

1 Laboratoire associe au C.N.R.S. 



of the Newton constant can modify the intrinsic geometry of the manifold. A phase 
transition is observed, separating two phases: in one the manifolds are crumpled, in 
the other they are elongated and resemble branched polymers. The transition was 
first discovered at d = 3, and is there of first order jjj. The latent heat is large, 
in agreement with the generic picture suggested by mean field arguments ||. The 
transition found at d = 4 was initially thought to be continuous ||, but a more 
recent study @ has furnished the evidence that it is actually of first order too, but 
with a very small latent heat. This finding has been confirmed by further work ||. 
The smallness of the latent heat remains, however, a mystery. On the other hand, 
the dynamics of the transition has been to a large extent elucidated []|, 10, pd[ |. 



In Ref. [12] it has been emphasized that the crumpled — > elongated transition 
occurs in all dimensions, but at d = 2 this happens only when the central charge c 
of matter fields is large enough. By lowering c below unity one enters the Liouville 
phase, where a sensible continuum limit can be defined. There exists a heuristic 
argument [|D|] explaining that: when, increasing c, one crosses the magic value 
c = 1 there occurs a condensation of metric singularities ("spikes") and a collapse 
of the Euclidean space-time. The continuum formalism and the known dynamics of 
the conformal factor (Liouville action) can be used to calculate the free energy of 
spikes. The authors of Ref. |l2| have suggested that a similar phenomenon might 
occur at d = 4. Using the effective action for the conformal factor calculated in 
Ref. [14], the free energy/spike isf]: 



1411 + N s + 11N F + 62AV - 28 
360 



In?. (1) 



The logarithmic factor on the right-hand side can be made arbitrarily large by taking 
a sufficiently spiky spike. The symbols Ns,fv refer to the number of massless scalar, 
fermion and vector fields contributing to the trace anomaly. The number 1411 comes 
from a one loop calculation of the contribution of transverse gravitons and 28 from 
ghosts and from quantum fluctuations of the conformal factor itself. This suggests 
that in the absence of matter fields the theory has no sensible vacuum and that this 
instability is lifted when the number of matter fields is large enough. As emphasized 
in Ref. []14| , it is remarkable that the sign of the contributions of matter and ghost 
fields is in 4d opposite to that found in 2d. Thus in 4d one stabilizes the theory by 
adding matter fields, while in 2d the contrary happens. This idea has been further 
expanded in Ref. |Tq] , 

In two dimensions the effective action behind the estimate of the contribution 
of spikes can be derived essentially rigorously. It should be emphasized that this 
is not the case in four dimensions, because it is not clear how to treat the trans- 
verse part of the metric. Nevertheless, we consider the conjecture discussed above 
sufficiently interesting to investigate the influence of matter on the geometry in 4d 
simplicial quantum gravity, although earlier investigations have not observed any 
non-trivial effects of the matter sector fl6[] . 

Since the right-hand side of Eq. (|l|) depends most strongly on the number of 
vector fields, we consider in this paper only continuous vector gauge fields, which 
have not been studied previously in 4d simplicial gravity. The model is defined in 
Section 2. In Section 3 we discuss the strong coupling expansion of the partition 
function for simplicial gravity interacting with any number of gauge vector fields. 
Our analysis of this series indicates that polymerization is suppressed when Ny > 2. 

2 The relevant spike configuration is <j) — ln{l + (? / \_(x — xo) 2 + a 2 ] }. The metric is propor- 
tional to and a <C p < const x a. 
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We check this point in Section 4, where the results of our Monte Carlo simulations 
are presented. Our speculations and conclusions form the content of Section 5. 



2 The model 

The action is a sum of two parts. The first is the Einstein-Hilbert action, which for 
a 4d simplicial manifold reads: 

S G = -k 2 N 2 + k a N a , (2) 

where denotes the number of /c-simplexes. The second part is 

Sm = E °(*ote) + Mhc) + A(l ca )f , (3) 

where A{l a b) is a U(l) gauge field living on link Z Q & and A{lab) = — A(Ji, a ). The sum 
extends over all triangles t a h c of the random lattice and o{t a bc) denotes the order of 
the triangle i Q & c , i.e. the number of simplexes sharing this triangle. Since we adopt 
a non-compact version of the theory, with a Gaussian action, there is no need to 
introduce a coupling in front of the sum on the right-hand side above. 

We work in a pseudo-canonical ensemble of (spherical) manifolds, with almost 
fixed N±. The model is defined by the partition function: 



AT 4 ) = J2 W ( T ) f' II dA V) e ~ SG ~ SM ~^ Na ~ (4) 

■' 1<ET 



The sum is over all distinct triangulations T and W(T) is the symmetry factor taking 
care of equivalent re-labelings of vertexes. The prime indicates that the zero modes 
of the gauge field are not integrated. As is well known, the volume conserving local 
move is not ergodic, hence we must allow the volume to fluctuate. The quadratic 
potential term added to the action ensures, for an appropriate choice of 5, that these 
fluctuations are small. The parameter k a is adjusted so that N A fluctuates around 
a mean volume N A . 

The Monte Carlo algorithm is constructed following one of the standard 



recipes of grand-canonical simulations [17]. A geometry move is done in two steps: 



first one integrates out the changing matter degrees of freedom, then the move is 
accepted/rejected based on a Metropolis test. If it is accepted, new matter fields 
are generated from a heatbath distribution, if new links are created, or disregarded 
if links have been removed. A geometry sweep consists of A4 attempted (randomly 
chosen) geometry moves. Matter fields are updated using heatbath and overrelax- 
ation algorithms. In a heatbath update the field is generated from the Gaussian 
distribution appearing in the partition function. In an overrelaxation update the 
field is given its "image" value: A(l) — * —A(l) + 2A, where A is the average over 
neighboring fields. A heatbath, or overrelaxation, sweep consists of N\ attempted 
updates. The relative frequency of geometry and matter sweeps is chosen so as to 
minimize both the auto-correlation time and the CPU demand; we use typically one 
heatbath sweep followed by one overrelaxation sweep, after each geometry sweep. 
Measurements are done when N4 = N4, and are separated by a constant number 
of "passes" through ./V4 = JV4. This number is chosen so that the system makes 
approximately 10 geometry sweeps between two successive measures. 
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Table 1: The number of different graphs, N g , for a fixed volume N4 and fixed number 
of vertexes iV , and the corresponding weights Wn v (N 4 , Nq). This is shown both for pure 
gravity, and one and three vector fields coupled to gravity (Ny = 0, 1, and 3). All weights 
are normalized with the value at N4 = 6. 
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3 The strong coupling expansion 

In parallel to the development of the Monte Carlo code, we have calculated succes- 
sive terms of the strong coupling expansion of the partition function Eq. ([|). This 
calculation is split into two distinct parts: First, for a given triangulation, one has 
to find the symmetry factor: W(T) = l/(number of equivalent re-labelings of T) . 
Second, one has to calculate the determinant A, resulting from the Gaussian inte- 
gration over one species of gauge fields (with unobservable zero modes kept fixed). 

We found the weights up to N4 = 18 by inspection, in order to test our codes 
and to get some insight into the problem. For larger manifolds we used pure gravity 
Monte Carlo simulations to identify the distinct triangulations and to determine 
the corresponding symmetry factors numerically, counting the relative frequency of 
the triangulations. The symmetry factors are also calculated explicitly, by going 
through all permutations of vertex labelings; comparing this to the Monte Carlo 
results serves as a consistency check on the identification. Finally, the determinant 
A is calculated using Maple. We have pushed the calculation up to N4 = 30; the 
different contributions to the partition function, Eq. (Q), are shown in Table 1 for 
zero, one and three vector fields. 

The big advantage of the strong coupling approach is that once the weights 
and the determinants have been calculated, it does not cost any effort to change K2 
and/or the number of matter fields. The series has been analyzed using the ratio 
method to extract the string susceptibility exponent 7 (cf. Ref. [0]), assuming the 
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Figure 1: Variations of 7 with K2, for 0, 1, 3, and 5 vector fields coupled to gravity. These 
values are obtains using the ratio method to analyze the strong coupling series, including 
terms corresponding to N4 = 6, 10, 14, 18, 22, 26, and 30. 



asymptotic behavior of the partition function: 

Z{N A ) ~ e K 4^4 i y7-3_ ( 5) 

Actually, only every second term of the series can be used for this purpose. One 
finds an oscillation in the behavior of the coefficients; a closer look reveals that 
the complete series is a sum of two. If one estimates the critical coupling using 
these series separately, the two estimates approach the correct value from below 
and from above, respectively. For 7 we obtained more reliable results using the 
terms corresponding to N4 = 6, 10, 14, . . .0 

Our results can be summarized as follows: 

(i) For fixed (large) K2 the susceptibility exponent 7 becomes negative when the 
number Ny of vector fields increases. 

(m) For fixed Ny and varying K2, the exponent 7 has a characteristic behavior: it 
is large negative for small enough K2, while it tends to a constant value for 
large K2- This limiting value expected, for pure gravity, and becomes 

negative for Ny > 2. For intermediate K2 the estimate of the exponent is 
highly unstable, which is presumably a reflection of the existence of a phase 
transition. This is illustrated in Fig. (|l]). For K2 > k§! 011 the other hand, the 
results are very stable, both with respect to variations in K2 and number of 
terms included in the analysis. 

These results suggest that for Ny > 2 the branched polymer phase disappears, 
and is replaced by a new phase with a negative 7. This is similar to what happens in 
two dimensions, except in four dimensions 7 becomes more negative as the number 
of vector fields is increased. All this is in qualitative agreement with the expectations 
described in the Introduction. As these results have been obtained using relatively 

3 A related regularity is observed in the distribution of the volumes of minimal neck baby uni- 
verses (minbu's) on the manifold. It has two branches, each corresponding to every second minbu 
size, which merge together at large minbu volumes. Also in this case, the series N4 = 6, 10, 14, . . . 
is less affected by finite size effects. 



5 



small triangulations, we have also carried out a series of Monte Carlo simulations 
to substantiate this picture; those results are presented in the next section. 

4 Monte Carlo simulations 

To understand better the phase structure of the model, as the number of vector 
fields Ny is increased, we have performed Monte Carlo simulations using one and 
three vector fields, mostly on lattices with no more than 16K simplexes. As these 
are rather modest lattice volumes, we present our results with the appropriate reser- 
vations; it is well known that for pure gravity, simulations with such small volumes 
can give misleading informations, especially about the nature of the phase transi- 
tion. We nevertheless believe that the results of our simulations, which agree well 
with the strong coupling expansion, present the correct qualitative picture. 

As we want to compare our results to the corresponding simulations of pure 
gravity, we summarize what is observed in that case: There is a first-order transition 

0, ^] separating a small K2 (strong coupling) crumpled phase from a large K2 (weak 
coupling) branched polymer phase; k\ ~ 1.29. The crumpled phase is characterized 
by two singular vertexes, connected to an extensive fraction of the total volume, 

1. e. their local volume grows like A4. Those vertexes are joined together by a sub- 
singular link; its local volume grows like A 4 . This singular structure is not present 
in the branched polymer phase, it dissolves at the phase transition. In simulations on 
small lattices one actually observes that the transition occurs in two steps II; 
first the sub-singular link is dissolved, later, at weaker coupling, the singular vertexes 
disappear. Those two sub-transitions, however, seem to merge as the volume is 
increased. 

This picture does not change much if one adds one vector field to the model. 
We still observe two peaks in the node susceptibility, xo = ((Aq) — (Ao) 2 ) / A4, pre- 
sumably associated with the two sub-transitions discussed above. Those transition 
separate the usual crumpled and branched polymer phases. The second transition, 
at larger K2, is more pronounced, and for the largest volume (32K) we observe a 
clear signal of a first-order transition, i.e. fluctuations between two states in the 
timeseries of the energy. The corresponding latent heat appears to be smaller than 
in the case of pure gravity — maybe an indication that the transition becomes softer 
as matter is added. 

A dramatic change occurs, however, when we introduce three vector fields to 
the model. We only observe a single peak in the node susceptibility, at K2 ~ 2. Its 
height increases with the volume; see Fig. (||a). Presently, our data are not good 
enough to decide whether this is a real phase transition and, in the affirmative, what 
is the nature of the transition and the corresponding critical coupling. Simulations 
with higher statistics and at larger volumes are needed for that. We also observe 
a change in the geometry of the manifolds coinciding roughly with the peak in xo'- 
at strong couplings the system is in the crumpled phase, one has two well identified 
singular vertexes of almost the same order, separated by a large gap from the orders 
of other vertexes. As K2 increases the orders of these vertexes merge with the rest 
of the vertex order distribution. This is shown in Fig. (||6), where is plotted the 
difference between the orders of the first and the third most singular vertex: po~P2- 
Combined, these results indicate that there still is a phase transition in this model. 

What about the respective phases? At small K2 there is the usual crumpled 
phase, but the weak coupling phase is no longer that of branched polymers. As 
mentioned above, it does not contain a singular vertex. This is best demonstrated by 
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Figure 2: (a) The node susceptibility xo, f° r three vector fields coupled to gravity, vs. 
K2- For 7V4 = 16000 the statistics is not sufficient for a reliable estimate of xo- (b) The 
corresponding change in the difference between the orders of the first and of the third most 
singular vertex, po and P2, normalized with the total volume N4. 



the scaling of the largest vertex order, which we show in Fig. (|3|) for K2 = 4.5. From 
the slope we see that the local volume grows like iV 4 3/4 , i.e. it is a vanishing fraction 
of the total volume. We have also measured the string susceptibility exponent 7 in 
this phase, extracted from the minbu distribution |]20|| . This we did for two values 
of K2 (4.5 and 6.0) and at different volumes (see Table 2). In all cases we get a 



Table 2: Measured values of the string susceptibility exponent 7 in the weak coupling (large 
K2) phase, for Ny = 3. 
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Figure 3: The scaling of both the largest vertex order, po, and of the average distance 
between two simplexes, ry, in the new weak coupling phase, for Ny = 3 and at k,i = 4.5. 
The lines are the best linear fits to the data. 



consistent negative value, 7 ~ —0.2, in reasonable agreement with the prediction 
from the strong coupling expansion, which is 7 ~ —0.38. The fits to the minbu 
distributions are very good, which confirms the validity of the assumed asymptotic 
behavior Eq. (||). In contrast, it is not possible to extract any reliable value of 7 in 
the crumpled phase, nor close to the phase transition. Finally, we have measured 
the Hausdorff or fractal dimension dn, using the scaling of the average distance 



between the simplexes [21 



<%>iV 4 = \X>^)\ ~ N l J d « , (6) 

were n(r) counts the number of simplexes at a geodesic distance r from a marked 
simplex. This is included in Fig. (||); a fit to Eq. (||) gives du = 3.97(15). 

One apparent pathology of this new phase, is that the node number is very 
close to its upper kinematic bound, i.e. (Nq)/N4 ~ 0.25. What this implies for the 
nature of the phase is not yet clear to us. 



5 Summary and discussion 

The most important result of this work is the discovery that by introducing matter 
fields one can prevent the standard collapse of 4d random manifolds into branched 
polymers. Contrary to all earlier investigations |TJ|, we find a strong back-reaction 
of matter on geometry. This opens new research possibilities although, of course, 
one still faces many uncertainties and pathologies. 

As already mentioned, the (discontinuous) phase transition observed with 
Ny = 1 seems softer than that in pure gravity. Eq. (|l|) indicates that a weighted 
sum of Ns,f,v controls the behavior of the theory; it is not excluded that, for some 
appropriate "combination" of matter fields, the transition might become continuous. 
This would be very exciting: one would have a relation between the matter content 
of the theory and its very existence. 

The nature of the new weak coupling phase, which we have observed, for 
Ny = 3, is also very intriguing. Our results suggest that this phase is characterized 
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by a non-trivial negative 7 and that it has fractal dimension four, the same as the 
flat space. Yet it has some pathologies. Although its largest vertex order scales sub- 
linearly with the volume, it grows much faster than, for example, the largest vertex 
order in pure 2d gravity (where the growth is logarithmic). In addition, we observe 
that the number of nodes of the manifolds is very close to its kinematic bound. In 
spite of this, it is not excluded that a non-trivial continuum limit can be taken in 
this phase, i.e. that the whole phase is critical, as a negative susceptibility exponent 
would suggest. That would be analogous to what happens in two dimensions, but it 
is not clear what would be the nature of such a continuum theory, obtained without 
tuning the Newton's constant to a critical value. 
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